Shaking-induced motility in suspensions of soft active particles 
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We investigate theoretically the collective dynamics of soft active particles living in a viscous fluid. 
We focus on a minimal model for active but non-motile particles consisting of N > 1 elastic dimers 
deformed by active stresses and interacting hydrodynamically. We first derive a set of effective 
equations of motion for the positions of the particles. We then exploit these equations in two 
experimentally-relevant cases: uncorrelated random internal stresses, and uniform monochromatic 
external shaking. In both cases, we show that small groups of intrinsically non-motile particles 
can display non-trivial modes of locomotion resulting from the hydrodynamic correlations between 
the particle-conformation fluctuations, fn addition, we demonstrate that a coherent shaking yields 
spatial ordering in suspension of soft particles interacting solely through the fluid. 
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I. INTRODUCTION 

Active suspensions generically refer to an ensemble 
of particles which continuously convert energy supplied 
by internal or external sources to deform and/or propel 
themselves in a fluid. Classical examples include sus- 
pensions of biological swimmers such as bacteria, algae 
or spermatozoa [1-4], microscopic robotic devices [5-7], 
and field-responsive particles which change their shape 
in response to external applied fields (such as electric, 
magnetic or temperature fields...) [8-10]. 

Recently, a number of theoretical studies have been de- 
voted to understand the large scale properties of this im- 
portant class of non-equilibrium systems [11-16]. Most 
of these coarse-grained models describe the large scales 
features of the active media, and ignore the specifics of 
the deformations or propulsion mechanisms at the sin- 
gle particle level. For instance, to model a suspension of 
swimming bacteria, the deformation pattern of the bacte- 
rias body is not explicitly described, but instead the cells 
are modeled as constant, local, and anisotropic stress dis- 
tributions acting on the surrounding fluid [11, 17]. The 
argument used to discard the details of the deformation 
dynamics is that these hydrodynamic models should be 
valid for time scales larger than the typical time-scale of 
the periodic microscopic actuation. 

However, this implicit time averaging step is not triv- 
ial, and it may in fact hide some subtle collective ef- 
fects [18-21]. Indeed, a priori, for a given active internal 
stress, the instantaneous rate of deformation of a soft 
active particle depends on the conformation, the orien- 
tation, and the position of all the other particles due to 
long-ranged hydrodynamic interactions. This leads to a 
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natural question to ask: Can the hydrodynamic coupling 
between the unsteady conformation and the particles po- 
sitions yield unexpected phase behavior and new modes 
of locomotion in suspensions of soft active particles? 

Here, we investigate a simple class of active suspensions 
comprising actively-deformed, but not-self-propelled, soft 
particles which interact only through the surrounding 
fluid. We start by rigorously deriving the effective equa- 
tions of motions for the center of the particles in the 
dilute limit. After discussing the case of deformation- 
controlled dynamics, we then consider the coupled con- 
formational and translational dynamics of the suspen- 
sion. We focus on two opposite, but experimentally rele- 
vant, situations. First, we assume that the particles are 
subject to stochastic but uncorrelated internal actuation. 
This is probably the simplest description to account for 
the shape changes induced by conformational-transitions 
in nano-machines [22] , and the shape fluctuations of non- 
motile living cells [23]. In that case, the center of a 
group of identical soft active-particles only experiences 
correlated random motion, but suspensions of heteroge- 
neous particles can display cooperative directed motion. 
In addition, finite-size suspensions evaporate or condense 
depending on their spatial distribution. Next, we con- 
sider the case where the particles are subject to a deter- 
ministic and uniform shaking at a prescribed frequency. 
This type of forcing is relevant to model the actuation 
of soft colloids, or more generally field-responsive parti- 
cles, by external fields. In that case, although isolated 
particles generate no net flow, hydrodynamic coupling 
induces nontrivial collective swimming modes for either 
identical or non-identical particles, as well long-range ef- 
fective interactions. Unbounded suspensions including 
a finite number of identical particles always evaporate, 
while polydisperse suspensions can also condense depend- 
ing on their relative spatial distribution. Remarkably, in 
the thermodynamic limit, systems made of identical par- 
ticles display a shaking-induced crystallization. 
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FIG. 1: Minimal model of soft active particle suspensions: 
N elastic dimers moving in the (x,y) plane. Top: 3D view. 
Bottom left: side view. Bottom right: top view. 



The paper is organized as follows. We introduce our 
minimal model of soft active particle suspension in sec- 
tion II. The general derivation of the effective equations 
of motion for the noise- or time- averaged positions of 
the N particles for the two opposite shaking sources (co- 
herent and incoherent) is then offered in section III. The 
case of incoherent (stochastic) forcing is studied in detail 
in section IV, while the case of coherent (deterministic) 
forcing is addressed in section V. We finish by some con- 
cluding remarks, and some suggestions for future work. 



II. MODEL 

Our minimal model is illustrated in Fig. 1. We con- 
sider TV dimers (what we will refer to as "particles"), 
each of them being made of two spheres of radius a con- 
nected by frictionless harmonic springs. We note ki and 
and £i(t) = £ + 6£i(t) the stiffness and the length of 
the i th dimer. All the particles have the same shape 
but can have different stiffnesses. The dimers live in a 
three-dimensional viscous fluid (Stokes flow), but are as- 
sumed for simplicity to remain in the two-dimensional 
(x, y) plane, and oriented along the z-axis, see Fig. 1. 
In that case, the dynamics of the system is described by 
the elongation £j and of the position of each of the N 
dimers. Each dimer in the suspension is force-free, but 
can expand and contract due to internal stresses modeled 
here by two equal and opposite forces acting on each of 
the two spheres. These forces of amplitude ff ct {t) are 
assumed to be oriented along the dimer axis, and the 
particles are thus deformed only in the direction normal 
to their motion. The linearity of the stokes equation im- 
plies that the velocities and the elongation rates of the 
dimers are linearly related to the total internal forces, 
which are the sum of an elastic (restoring) and of an ac- 



tive components, fi(t) — ff(t) + /f ct (i), and thus 
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In Eqs. (1) and (2), fj refers the z-magnitude of the 
force acting on the top dimer in Fig. 1 (bottom left). 
In principle all the mobility coefficients m,j and in 
Eqs. (1) and (2) depend both on the elongation of all 
the particles, as well as on all the inter-particles dis- 
tances, Tij{t) = dij + Sdji(t). This naturally makes the 
above equations highly nonlinear, and in general impossi- 
ble to solve by hand. In order to obtain analytical results 
we restrict our study to the limit of dilute suspensions, 
dij/£i 3> 1, and wide dimers, £/a 3> 1. In that limit, the 
2N spheres composing the TV dimers can be replaced by 
point forces (Stokeslets). The velocity of the n th sphere 
is thus related to the force acting on the m th sphere by 
<9 t R„ = Ysm G ( R n ~ Rm) • F m (t), where R„ is the posi- 
tion of sphere n, and G = (l + RR/i? 2 ) /(8nriR) is the 
Oseen kernel [24]. Within this approximation, the mo- 
bility coefficients defined in Eqs.(l) and (2) can be com- 
puted analytically at the necessary order in in a/£ and 
£/r. Setting G(r) = l/(47rryr), and G(0) = 1/(6^), it 
is straightforward to get 

li u = 2G(0) = /u , 
(Hj = -j^iG"(r y ), i ^ j 



m, 



m, 



(3) 
(4) 
(5) 
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where is the unit vector pointing from the center of 
mass of dimer i to that of dimer j, and with higher-order 
terms which are negligible in what follows. Importantly, 
we note that fj, is symmetric and that m is not antisym- 
metric. At higher order in the separation distance be- 
tween the dimers, only corrections with even (resp. odd) 
derivative in G add up to the /iy (resp. to the m^). 
Finally, we note from Eqs. (1) and (5) that the position 
of an isolated dimer does not vary in time: The particles 
are therefore active, but are not self-propelled. This re- 
sult is a direct consequence of PurcelPs scallop theorem 
stating that body deforming with time-reversible strokes 
cannot display any average motion in a Stokes flow [25] . 



III. EFFECTIVE DYNAMICS 

In this section, we derive the coarse-grained equations 
of motion for the dimers. Specifically, we explain how 
to adequately perform averages for Eq. (1). All along 
this section, we use the symbol (...) to denote averages. 
Anticipating the two types of actuation of interest, these 
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averages can represent two different things, fn the case 
of a stochastic forcing, the averages are made over a large 
number of noise realizations (section IV) . In the case of 
a deterministic AC forcing, they refer to time-averages 
over the actuation period (section V). The algebra in 
both cases is identical, so we present it below without 
specifying which of the situation we are referring to. The 
final equation we obtain is Eq. (16), and the reader un- 
interested in technical details can skip what precedes it. 



and By. At lowest order in S£/d, Sd/d and a/l we have 
A i:j = my ^^(SdjkSik) + my ^ ^lL(S£jS£ k ), 



and 



(11) 

(12) 



Derivation of the coarse-grained equations of 
motion 



We restrict ourselves to the limit of small deforma- 
tions, S£i/£i <C 1 and small displacements Sdij/dij <C 1. 
Expanding the mobility coefficients, my at first order in 
S£ and averaging Eq. (1) yields 



(d t ri) 
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m ij(fj) 



dd i:j 



(Sdijfj) 
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n=i,j 
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where Taylor expansions have been used, and therefore 
all the my terms in Eq. (7) are evaluated at £i and dij. 
Importantly, in Eq. (7) and throughout the paper, we 
will not be using summation over repeated indices, and 
therefore sums will be explicitly written down. 

We are now going to show that in the dilute limit, the 



where we have used (5£ n ) = and (S£iS£i) — 0. Since 
we have dy = (rj — ri) 2 , we can use Eqs. (7) and (8) to 
write Sdij = eij ■ (Srj—Sri) at first order in the elongation 
variables, and therefore get 
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Using Eq. (13) in Eqs. (11) and (12), and recognizing 
again that (S£iS£i} = 0, we see that the coefficients in 
front of averages of the form (6£j5£k) in both Ay and By 
are of the order of m times an off-diagonal term of the 
matrix and therefore decay in space as 0(G'G") ~ 
1/r 5 for an unbounded fluid. 
We now turn to the Cy term. At lowest order we have 



Cy — /i 



(14) 
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^G"'(dy)<*W->e«. (15) 



Consequently, in the Cy term, the coefficient in front 
of the conformation averages are on the order of G'" , 



right hand side of Eq. (7) is dominated by the Cy terms. which decays in space as ~ 1/r 4 in an unbounded fluid. 



To identify this dominant contribution, we expand the 
forces from Eq. (2) at second order in the 5£i and Sdij 
and obtain 



Si 




t)ft 



lS£ k + -^Sd jk S£ k +J2^j L S£ n S£ k ] ,(8) 
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n=j,k 



8£n 



Since G'" 3> G'G" , we see therefore that, in the dilute 
limit, the Cy term in Eq. (7) dominates the dynamics of 
the time-averaged positions of the dimers. Consequently, 
Eq. (7) reduces to the final coarse-grained equation 



d t {vi) 
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Y,G"'{dij){S£i5£j)eij, (16) 



where At„„ refers to the ran component of the matrix fx~ l 
inverse of /i. We then use this expansion to compute the 
Aij, Bij and Cij terms in Eq. (7). 

To do so, we fist need to compute the inverse of the 
mobility matrix /x. In the dilute limit, we get from 
Eqs. (3) and (4) that fm = fi ^> % for all i and 
j =/= i. As a result, we can expand [i at first order in 
the Hij/no ~ £ 2 G"(dij)/G(0), and invert the matrix to 
obtain at the same order 
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Let us now first consider the order of the terms A, 



which is our first main result. 



B. Deformation-controlled scenario 

For a system composed of more than one particle, 
Eq. (16) shows that the dimers positions will evolve in 
time if their conformations have nontrivial correlations. 
So, a collection of non-swimmers can be made to swim. 
Physically, this cooperative motion arises due to the cou- 
pling between the location of each particle and the instan- 
taneous dipolar flow field induced by the deformations of 
the surrounding particles. 

In the case where the conformation kinematics is pre- 
scribed, then Eq. (16) fully prescribes the suspension dy- 
namics. For instance, let us consider the simple case of 
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two identical dimers with sinusoidal length modulations: 
£i(t) = £ + 5£cos(ujt + 4>i). It is then straightforward to 
derive from Eq. (16) the overall average swimming ve- 
locity, v = i(<9 t (ri + r 2 )), and overall relative velocity, 
Av = (9t(r 2 — ri), and we find 



27 (£5t 



rf 2 



L)ei2, 



Av 



0. 



(17) 
(18) 



The two dimers remain at a constant distance on average. 
If <p2 — (f>i — or 7r, then the deformations of the two- 
particle system are reciprocal, and average swimming is 
prohibited by the scallop theorem. In all other cases, 
the two particles swim collectively. There is therefore no 
many-scallop theorem [20, 21]. 

Next, and more formally, it is important to recog- 
nize that the effective dynamics of N soft active par- 
ticles cannot be mapped on the relaxation dynamics of 
an ensemble of passive particles interacting through a 
long-ranged pair potential. This is because the confor- 
mation correlation terms in Eq. (16) has the proper- 
ties of a pseudo-scalar. Indeed, the length fluctuations 
are either periodic or statistically stationary, which im- 
plies that {5li{t)5£j{t)) = -{5£j{t)5£i{t)). In addition, 



G"'{d 



VG"(dij) derives from a scalar quantity, 
this implies that the full term in Eq. (16) cannot derive 
from a scalar potential, and our conclusion is therefore 
true for all type of actuation applied to the soft particles. 

In the sections below we focus on two different, and 
more complex, types of actuation. First we consider the 
case of particles shaken by internal stochastic forces hav- 
ing no-spatial and no-temporal correlations (incoherent 
shaking, section IV). Next we study the case of a uniform 
sinusoidal external force acting on all the soft particles 
(coherent shaking, section V). In both cases, we will solve 
Eq. (2) whithin a linear response approximation to com- 
pute the conformation correlations, (S£iS£j), and Eq. (16) 
will then be used to derive the effective dynamics of the 
suspension. 



IV. INCOHERENT SHAKING 

A. Effective dynamics 

In this first scenario, we assume that the active com- 
ponent of the force on each sphere in a dimer is a random 
variable with zero mean and fluctuations given by 



(.ff ct (i)> = 0, 

(/rw/n*')) = frs(t 



t')Sij, 



(19) 
(20) 



where / is the typical force amplitude and r _1 the high- 
est frequency of the force power spectrum. Since / act is 
not spatially correlated whereas the mobility matrix fi 



is non-local (see Eq. 4), the fluctuation-dissipation rela- 
tion cannot for this noise term, and therefore / act does 
not represent thermal noise. Instead, the term / act mod- 
els active and fluctuating small-scale forcing, for example 
the stochastic actuation by molecular motors. 

To derive the conformation correlations we separate ex- 
plicitly the active component of the force from its elastic 
component in Eq. (2) 



dt.S£i 
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Hijkj5£j + Hijfj (t), 



(21) 



where we evaluate the fiij(£i 7 £j,rij) term at (£,£,dij). 
Within this linear response framework, we Fourier trans- 
form Eq. (21) to obtain 



S£i(u)) = ^Rij^fjiu), 



(22) 



where R(w) = [iwl + /xk(w)] -1 fi{uj) is the kernel given 
by Eq. (21), and the matrix k is defined by fcjj = hSij. 

By writing the Fourier-transform of Eq. (20) as 
(/f ct (w)/j Ct (u/)> = 2tt/ 2 t(5(w + uo')8 lJ1 we can then use 
Eq. (22) to obtain the correlation as given by 

(5Ut)5£ 3 (t')) = (23) 

k 

We next calculate the coefficients of the response matrix 
R for dilute suspensions. In this limit, the off-diagonal 
terms of fi are subdominant, and after some straightfor- 
ward algebra, and at the lowest order in the we 
get that the Rij terms reduce to 



Mo 



i + j, 



(24) 
(25) 



where we have defined the relaxation frequency of the 
dimers by u>i = /io^i- The sum over k in Eq. (23) can 
next be approximated by the terms at the leading order 
in 1/ dij . The dominant contribution come from the two 
terms in the form of a product of diagonal terms in R 
times an off-diagonal term, i.e. 



Y,[R ik (uj)R jk (-(j)] = 

k 

Rui^Rjii-uj) + Rij^Rjji-u) + h.o.t. 



(26) 



where the higher-order terms decay faster in space. Using 
Eq. (24) and (25), and after integration of Eq. (23) over 
ui we get the correlation as given by 



<<%(*)«,•(*)> = J(A)mo/%-^-^' 

Z UJi -f- U)j 



(27) 



Finally, using the definiton of the Eq. (4), and the 
above expression for the correlation in the equation of 
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motion Eqs. (16), we obtain the final expression for the 
effective equations of motions of the N soft dimers as 

dtin) = £ ^^G"{d ij )G"'{d ij )e ij . (28) 
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B. General comments 

The general equation describing the average motion of 
each soft active particle is Eq. (28). Interestingly, we 
first note that this final formula does not involve the 
conformational self-mobility coefficients of the individ- 
ual particles, /ip, emphasizing the role of hydrodynamic 
interactions. In particular, the same final formula would 
be obtained for dimers of different shapes and sizes. The 
fast spatial decay of the velocities in Eq. (29) with the 
inter-particle distance should also be stressed. In an un- 
bounded fluid, G ~ 1/d and therefore G"G"' ~ 1/d 7 . 
This is to be contrasted with the G'" ~ 1/d 4 power 
law predicted when the deformations are imposed (see 
Eq. 16). This difference arises from the 1/d 3 decay of the 
elongation-elongation correlation function in a force con- 
trolled set-up, whereas the correlation is by construction 
independent of the distance when the deformations are 
imposed. 



C. Small groups of soft active particles: N = 2 

The right-hand side of the equations of motion, 
Eq. (28), is the sum of the N pairwise hydrodynamic in- 
teractions between the dimers. Some additional insight 
on the collective dynamics can thus be gained by consid- 
ering first the case of a single pair of dimers. For N = 2, 
it is easy to see in that the average and relative speeds 
are given by 





v16tt) ' 
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+ UJ2/ 



/VI 4 

f] 2 d 7 



ei2, 



Av = 0. 



(29) 
(30) 



The two particles remain at a constant distance for all 
noise amplitude. In fact, if the two particles are identical 
(same relaxation time) they do not move at all in aver- 
age (v = 0) . If however they posses different relaxation 
times, the particles move together with the same speed 
oriented along the ei2 direction, and with a sign that 
depends on the relative magnitude of their relaxation 
frequencies. This collective mode of locomotion arises 
despite the fact that (i) the two particles are non-motile 
when isolated and (ii) the stochastic shaking forces do 
not induce any averaged fluid velocity in the far field. 



D. Interpretation 

To understand intuitively the existence of a net lo- 
comotion in the presence of noise, we can think of 



the stochastic shaking forces as the incoherent super- 
position of sinusoidal forces with uncorrelated random 
phases. The response to two forces of the form ff ct — 
f J ducos[ujt + ipi(u))] results in out-of-phase oscillations 
of the two dimers lengths, with a phase difference hav- 
ing two contributions. The first contribution comes obvi- 
ously from the phase differences between the two shaking 
forces; the second intrinsic contribution arises from the 
difference between the two relaxation times. 

Let us first ignore the off-diagonal coupling of the con- 
formation mobility tensor fj^,j. In that case, the dimer i 
would respond solely to the force and would thus os- 
cillate with the phase u)t + ipi + fa , where fa is the phase 
shift due to the finite response time of the elastic spring 
coupled to the viscous fluid. In turn, the elongation of the 
first dimer and the elongation rate of the second dimer 
would be necessarily uncorrelated as the shaking forces 
have uncorrelated phases. More precisely, the average ve- 
locity given by Eq. (17) scales with sin(</>i +ipi — fa — "02 ), 
which averages to zero for any two relaxation times. 

If we now take into account the off-diagonal hydrody- 
namic coupling, we can physically interpret the appear- 
ance of locomotion. Indeed, as momentum propagates 
instantaneously in the fluid (there is no inertia), the elon- 
gation of, say, the first dimer, i — 1, has now two contri- 
butions. The first contribution is the local response to the 
local force /i, while the second results from the dipolar 
flow field which oscillate in phase with the dipolar force 
sources, and is induced by the force f 2 acting on the sec- 
ond dimer. The local contribution to the elongation has a 
phase given by ut+ipi + fa, whereas the non-local contri- 
bution to the deformation is oscillates with cut + fat + fa ■ 
It then follows that the local contribution to the elonga- 
tion of the first dimer correlates with the non-local con- 
tribution to the elongation rate of the second dimer if the 
two particles have different relaxation times. Indeed, this 
correlation scales now with sin(fa +ipi — fa — ipi) (and an 
amplitude decaying with d), which does not average to 
zero as the dimers are distinct. Note that such interpre- 
tation also holds for viscoelastic dimers which internally 
dissipate energy (results not shown). 



E. Small groups of soft active particles: N = 3 

Before addressing the dynamics of large N system, we 
briefly consider the effect of the presence of a third parti- 
cle (N = 3 case). For simplicity we assume the particles 
to be located along a line (x axis), and restrict the hydro- 
dynamic kernel G to nearest-neighbor interactions. By 
scaling times by ( 51 ^ / 2 t1 4 /?7 2 ) and numbering the 
particles along the positive x axis, Eq. (28) becomes 



u>i — U>2 



9 t Xl -- ; | „ 

+ o>2/ d[ 
d t x 2 = d t x 3 + d t x\. 



uj 2 - 0J3 
w 2 + W3 / dl 



,(31) 
(32) 
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If we denote d\ = d\ 2 = x 2 — X\ and d 2 = d 23 = %3 — x 2 , 
it is easy to see that the inter-particle distances obey 



d t di = [ ^ 2 ; ^ 3 ) -L, d t d 2 



lu 2 + u> 3 J dl 



^^)~ (33) 

Wl + LJ 2 J d\ 



Obviously the only interesting case is the one where the 
three intrinsic frequencies are all different. Without loss 
of generality we can assume lu\ < lu 3i and the only pa- 
rameter is the relative magnitude of lu 2 with respect to 
these two frequencies, so generically three different cases 
exist. If uj 2 < ^1 < ^3 then the suspension collapses 
(condenses) in finite time, and swims only on a finite 
length. If uji < uj 2 < u 3 then one side (here di) collapses 
in finite time as <~ (to — t) 1 / 8 and the other side (here, d 2 ) 
diverges (evaporates). If < u> 3 < u 2 the whole system 
evaporates as di ~ i 1 / 8 , leading to an average velocity 
on the order of v ~ i~ 7 / 8 , and swimming over an infinite 
distance Ax — J vdt <~ i 1 / 8 . 



F. Large groups of soft active particles: N 3> 1 

The general equation describing the average motion 
of each soft active particle is Eq. (28). The right-hand 
side of this equation is the sum of the N pairwise hy- 
drodynamic interactions between the particles. Since all 
two-body interactions vanish when the relaxations times 
are identical (u)i =u)j), a suspension of identical soft ac- 
tive particle does not display any net motion when driven 
by incoherent active stresses. Any spatial distribution of 
the particles would thus correspond to a stationary state 
(d t (ri) = 0). Nonetheless, the ensemble of soft active 
particles should display collective diffusion, the dynam- 
ics of which could be in principle obtained following the 
same framework as the one used above. 

Oppositely, suspensions of non-identical particles have 
in general a non-zero collective velocity. If an en- 
semble of such particles is randomly distributed in the 
plane, with the particles having a non-trivial (and well- 
behaved) distribution of relaxation time, it follows from 
the central limit theorem that the collective velocity 
V = ^Adtfi) /N has a random orientation, and a fi- 
nite amplitude that scales as V ~ l/y/~N. In addition, 
the direction taken by such a "swarm" of soft active par- 
ticle would be dictated by the initial anisotropy of the 
particle distribution. 



V. COHERENT SHAKING 

In this section, we consider the scenario where we have 
a coherent driving. The active component of the force 
is assumed to be applied identically on all dimcrs, and 
to vary sinusoidally in time with prescribed frequency u, 
i.e. ff ct {t) — fcos(u)t). This case is relevant to a wide 
range of experimental situations where a uniform exter- 
nal AC field is used to change the shape of responsive 



particles suspended in a viscous fluid. Examples include 
magnetic fields used to bend arrays of magnetic colloidal 
filaments [9] , light or temperature cycling used to change 
the shape of nematic elastomer colloids [10], or electric 
fields used to deform emulsions of dielectric droplets. 

Within the same linear response scheme as the one 
discussed in section IV, the deformations of the particles 
oscillate at the same frequency as the forcing, ui, and are 
given by 

5£i(t) = 8£iCos(ut + <j>i), (34) 
where the amplitude and phase are given by 



Mi = f 



arg 



(35) 



In order to compute the effective equation of motion, 
Eq. (16), we need to calculate the time-averages 
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(SliSij) = -u\5ti\\5£j\ sin(0 4 - 



(36) 



and the analysis differs between suspensions of identical 
or non-identical particles. 



A. Suspensions of non-identical particles 

We first consider the situation where the soft parti- 
cles have different relaxation times. In that case, and in 
the dilute limit, the response matrix R is dominated by 
its diagonal terms. Exploiting Eqs. (24) and (25) it is 
straightforward to get 



m\ = 



Mo/ 



and therefore 

nlf 2 uj 2 {u t - ujj) 



i = — arctan [ — ) , (37) 



31 2 (L0 2 +Luf)(Lu2+Lu])> 

so that the effective dynamics is fully described by 

dt{Ti) 



(38) 



31 Wo u 2 (uj-ujj) 



32 (a ;2 + w 2 )(w2 + w 2 ) E G/ "(^) e «- 



(39) 

The result of Eq. (39) generalizes to N particles the anal- 
ysis done in Rcf. [21] for TV = 2 dimcrs. We obtain an 
effective dynamics with a <~ 1/d 4 spatial decay. 



B. Suspensions of identical particles 

We now focus on the case of an ensemble of identical 
soft particles. We assume that the particles all have the 
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same stiffness, k, and therefore identical relaxation fre- 
quency, LOi = lu . The amplitudes in Eq. (36) are still 
dominated by the diagonal terms in R and we have 



|<%| = 



Mo/ 



(o;2+u; 2 )i/2 



(40) 



To get the phase differences, we have to look at the con- 
tribution of the off-diagonal terms in R, as the leading- 
order term for fa given by Eq. (37) predicts no phase 
difference for identical dimcrs. Exploiting the fact that 
we are in the dilute limit, and that therefore the phase 
differences are expected to be small (<j>i — (fij <C 1), we 
can use Eqs. (24) and (25) to explicitly obtain at leading 
order in 1/dij 



- 03 = 



1 UJLUq 
[If) LU 2 + LJq 



k^i k^j 



(41) 



The effective suspension dynamics is then obtained to be 
U 2 .f 2 lu 2 lu 



dtiri) = 



32 {uo 2 + w 2 ) 



2\2 



k^i kjtj 



eij, (42) 



where the {fi mn } are given by Eq. (4). 



C. General comments 

As a difference with the results of section IV, we see 
that in general, identical particles subject to AC forcing 
do display collective dynamics. Since only off-diagonal 
terms of fi appear in Eq. (42), the spatial decay of the 
velocity field scale generically as G"'G" <~ 1/d 7 in an un- 
bounded fluid. We also note that Eq. (42) actually rep- 
resents a three-particle interaction, with a non-pairwise 
additivity of the effective coupling. For a homogenous 
suspension with N = 2 we have dt{x\) = dt(x2) = 0, and 
no motion occurs on average, whereas N — 3 particles 
do move (see below). To provide with a physical inter- 
pretation of the results quantified by Eq. (42), we now 
proceed to analyze two cases in detail. 



be all the more valid that the system will be confined. Us- 
ing the explicit expressions for the mobility coefficients, 
Eqs. (3) and (4), we then obtain the form of the equations 
of motions for each dimer as 

dt(xi) = A(iv)G'" (di) [G"(d i+1 ) - G"(di-i)] 

+ A(u>)G'" (di-i) [G"(di) - G"(di-2)] , (43) 

where we have defined di = di.i+i and where the pref- 
actor A(w) is given by A(w) = -^£ 4 f 2 uj 2 ujo/(u 2 + lo 2 ) 2 . 
Importantly, Eq. (43) shows that changing the forcing 
frequency does not induce any qualitative changes in the 
dimers dynamics as it can be absorbed into a rescaling 
of time. We are thus left with a set of purely geometrical 
equations of motion. 

An obvious stationary state of Eq. (43) corresponds to 
a crystal order where all the distances are identical, di = 
d, To investigate the stability of this ordered stated, 
we apply a small perturbation to the particle separations, 
di = d + Sdi(t), and linearize the equations of motion to 
get 

d t Sdi = AH [G"'(df {Sd i+2 - 2Sdi + Sdi- 2 ). (44) 

To compute the eigenvalues of Eq. (44) we look for wave- 
like solutions of the form 5d n (t) = 5d n exp[z (ilt + Qnd)] 1 
which leads to the dispersion relation of the compression 
waves along the chain of active dimers 



iVL = -4A(w) [G"'(d)} sin 2 (Qrf), 



(45) 



with similar results obtained if the nearest-neighbor ap- 
proximation is relaxed. For all shaking frequencies, the 
compression modes are overdamped, the crystal state is 
stable. In the low Q (wavenumber) limit, the above 
dispersion relation is identical to the one of a chain of 
particles damped by a local drag force, and interact- 
ing through a harmonic potential. In other words the 
shaking induces an effective elasticity to this array of N 
non-interacting particles only dynamically coupled by the 
surrounding fluid flow. Note that such mapping on the 
relaxation dynamics of a passive system is correct only 
for small perturbations around the steady state, and is 
in general not appropriate to describe the motion of the 
dimers (see section IIIB). 



D. Shaking-induced crystallization 

We first consider the dynamics of infinite one- 
dimensional chains. We assume that the soft particles are 
all aligned along the x direction, and number them such 
that Xi < Xj if i < j. Such a setup may be achieved ex- 
perimentally by studying the dynamics of field-responsive 
particles confined inside a capillary tube or in a micro- 
channel. As the ii mn coefficients in Eq. (42) decay as 
1/d 3 , we further make a nearest-neighbor coupling ap- 
proximation and neglect any {i, i+n} hydrodynamic cou- 
pling with \n\ > 1. Experimentally this assumption will 



E. Small groups of soft active particles: 
Evaporation and spontaneous locomotion 

We now focus on finite-size one-dimensional chains 
(with N > 2). In that case, the suspension docs not 
posses a stationary state. This is best seen by writing 
down the evolution equation for the two end-particles, 
which are, under the nearest-neighbor approximation, 

d t ( Xl ) - K{w)G'"{d l )G"{d 2 ), (46) 
d t (x N ) = -K{uj)G"'{d N - 1 )G"{d N - 2 ). (47) 
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We see therefore that dt(x\) < as well as dt(xjv) > 
0, and therefore for any conformation of the chain, the 
velocity of the two end-particles points outwards, The 
finite chain of active particles is thus seen to evaporate. 
Such evaporation takes place in all finite suspensions, as 
confirmed by numerical simulations in two dimensions 
(not reproduced here). 

The evaporation dynamics can be analyzed in more 
detail in the case with N — 3. In that situation, using 
the exact expression for G(d) and rescaling time to absorb 
all irrelevant prefactors, the three equations of motion 
reduce to 

dt{xi) = -rfPf' dt{X3) = rfPf' (48) 
d t (x 2 ) = d t { Xl )+d t {x 3 ). (49) 

It is then straightforward to show that (d t di)/di = 
(dtd 2 )/d 2 and thus that the ratio between the two dis- 
tances, a = d 2 /di, is a constant of the motion. It then 
follows that the inter-particle distances grows with time 
as d x {t) ~ a- 1 /2 t V8 and d 2 (t) ~ a 1 / 2 * 1 / 8 . Hence the 
growth speed of the suspension varies as d t (di + d 2 ) ~ 
(a 1 / 2 + a -1 / 2 )* 1 / 8 . Let us now consider the average 
speed, v = ^(d t {xi) + d t (x 2 ) + d t (x 3 )), of the three 
dimers. It is easy to see that v = | (1 — a)d t di. For 
a symmetric conformation, di = d 2 , and the system does 
not swim on average. However, the hydrodynamic cou- 
pling between the dimers induces a spontaneous left-right 
symmetry breaking. Indeed, an infinitesimal perturba- 
tion to the ratio, say d 2 jd\ = 1 + e, cannot relax as 
d 2 jd\ is a constant of motion. This small swarm of ac- 
tive particles acquires therefore a spontaneous net mo- 
tion with a velocity decreasing with time as v ~ t~ 7 / 8 , 
implying that the center of the suspension can travel 
over an infinite distance, Ax ~ i 1 / 8 . Three identi- 
cal active soft particles swim therefore cooperatively. 
Note that the dilatation width, d t {{d\ + d 2 )), and the 
swimming speed, v, decrease with the same scaling, and 
d t ((di + d 2 ))/\v\ = 3(1 + a)/2\l - a\ > 1, which implies 
that the system evaporates faster than that it swims for 
any initial conformations. 

F. Large groups of soft active particles 

Within the approximations made in this paper, the col- 
lective mode locomotion discussed above in the case of 
N = 3 can only exist in the case of finite suspensions. In- 
deed, from Eq. (43) it is straightforward to show that for 
N particles, the average speed of the suspension, v, scales 
^ v ~ jf - up to 0{l/N) bound- 

ary terms. So, if the di are finites, we have the inequal- 
ity jf\d N ~ d 1 |/max(d l ) 8 < v < j^\d N - d x \/ min(d 4 ) 8 . 
As a result, for all spatial distributions of the elastic 
dimers, the collective locomotion speed if at most of order 
0(1/N). Whereas spontaneous motion takes place for 
small groups of soft active particles, the effect vanishes 



therefore in the thermodynamic limit. Note however that 
a non-zero collective locomotion could arise from nonlin- 
ear effects in the conformational dynamics of the parti- 
cles. This non-trivial extension lies beyond the scope of 
this paper. 



VI. CONCLUSION 

We now conclude by returning to the initial question 
that motivated this work: Does the hydrodynamic in- 
terplay between the shape fluctuations and the positions 
of soft active particles yield unanticipated collective be- 
havior? As we demonstrated in this paper, the answer 
is yes. On the basis of a minimal model for hydrody- 
namically coupled soft dimers, we showed that despite 
the fact that these particles are intrinsically non-motile, 
we obtained the following results for their collective dy- 
namics. Finite ensembles of non-identical soft particles 
experience collective swimming when driven by an uncor- 
rected stochastic internal forcing. The direction and the 
amplitude of the resulting collective swimming is set by 
the initial anisotropy of the particles positions. In addi- 
tion, these "swarms" condense or evaporate depending on 
the distribution of relaxation times in the suspension. In 
contrast, soft particles having identical relaxation times, 
can only display collective diffusion. When shaken by a 
monochromatic and homogeneous field, finite ensemble 
of identical soft active particles always display collective 
swimming and evaporation in the dilute limit for all ini- 
tial positions. In the large N limit, the suspension dis- 
plays shaking-induced crystallization. However, in this 
large TV limit, the collective swimming velocity vanishes 
algebraically due to the fast decay of the correlation of 
the shape oscillations (1/d 3 ). 

The cancelation of the collective swimming speed in 
the thermodynamic limit might result from the linear re- 
sponse approximation we used to compute the shape fluc- 
tuations of the particles. In the dual and simpler problem 
of collective pumping, it has been recently shown that an 
array of hydrodynamically-couplcd filaments driven with 
the same constant torque can be unstable, resulting in 
the propagation of metachronal waves, and fluid pump- 
ing, see e.g. [26, 27]. A similar rectification phenomena 
could in principe exist in suspensions of freely moving 
soft active particles. This intrinsically non-linear effect 
would likely result in a motility transition in suspensions 
of reversibly actuated particles. 
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